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Technical  Summary: 

Summary  of  Papers 

Goldberg,  Z.,  Rocke,  D.,  Schwietert,  C.,  Berglund,  S.,  Santana,  A.,  Jones,  A.,  Lehmann, 
J.,  Stem,  R.,  Lu,  R.,  Hartmann-Siantar,  C.  (2006)  Human  in  vivo  dose-response  to 
controlled,  low-dose  low  linear  energy  transfer  ionizing  radiation  exposure.  Clin  Cancer 
Res.,  12,  3723-3730. 

This  paper  is  the  first  to  describe  transcriptional  responses  to  low-dose  ionizing  radiation 
in  vivo  in  humans.  The  earlier  dosimetry  work  allowed  for  accurate  quantitation  of 
precisely  calibrated  doses  which  were  verified  by  TLD’s  and  MOSFETs.  Eight  men 
were  included  in  this  study,  a  biopsy  was  taken  prior  to  each  patient’s  first  radiation 
treatment  and  three  biopsies  were  taken  3  hours  post  exposure  at  locations  on  the 
abdomen  which  received  1 ,  1 0  or  1 00  cGy  doses.  Each  biopsy  was  processed  and  the 
RNA  was  interrogated  on  an  Affymetrix  HGU133  Plus  2.0  array.  The  data  was  analyzed 
using  the  method  of  Rocke  et  al  as  described  above,  which  determined  significant 
differential  expression  of  pre-selected  gene  groups  and  pathways.  The  results  identified 
seven  of  nineteen  gene  groups  and  five  of  seven  gene  pathways  which  were  up  or  down 
regulated.  The  gene  groups  up  regulated  included:  BCL6,  cytokines,  mitogen-activated 
protein  kinases,  and  zinc  finger  proteins.  The  gene  groups  down  regulated  included: 
keratins,  protein  disulfide  isomerase  and  SI 00.  The  five  pathways  all  displayed  up 
regulation  and  included:  transforming  growth  factor  p/cyclin/ubiquitin  pathway, 
stress/apoptosis,  inflammation,  growth  factor/insulin,  and  AKT/phosphoinositide-3- 
kinase.  The  gene  groups  and  pathways  found  in  this  study  have  known  functions  in 
apoptosis,  cell  cycle  arrest,  DNA  repair  and  inflammation.  This  indicates  that  the  skin 
responds  to  the  injury,  initiates  a  stress/inflammatory  response,  enhances  survival  and 
begins  DNA  damage  repair.  These  data  show  that  there  is  a  definable  response  in  three 
dimensional,  full  thickness,  in  vivo  tissue  to  LDIR. 


20080108188 


Berglund,  S.R.,  Rocke,  D.M.,  Dai,  J.,  Schweitert,  C.,  Santana,  A.,  Stem,  R.,  Lehman,  J., 
Hartmann  Siantar,  C.,  and  Goldberg,  Z.  (2007)  Transient  genome-wide  transcriptional 
response  to  low-dose  ionizing  radiation  in-vivo  in  humans,  International  Journal  of 
Radiation  Oncology,  Biology,  Physics,  in  press. 

This  is  our  second  published  paper  using  our  developed  methodology  and  is  the  first 
report  to  describe  the  temporal  transcriptional  changes  found  in  response  to  low  dose 
ionizing  radiation  in  vivo  in  humans.  A  different  cohort  of  men  was  used  in  this  study 
and  biopsies  were  taken  at  10  cGy  location  before  treatment  at  time  0  and  at  3,  8  and  24 
hours  post  exposure.  We  compared  the  0  and  24  hour  samples  to  the  3  and  8  hours 
samples  to  determine  transient  responses.  To  summarize,  we  assumed  that  responses 
seen  in  the  3  and  8  hour  biopsies  could  be  considered  transient  responses  if  they  returned 
to  the  baseline  seen  in  the  control  biopsy  by  the  24  hour  time  point.  Using  the  same  gene 
group  and  pathway  probe  sets,  we  found  that  9  gene  groups  showed  transient  expression. 
Groups  that  were  up-regulated  included:  zinc  finger  proteins,  keratins,  BMP  receptors, 
BAG,  and  cyclins.  Down-regulated  groups  included:  TNF,  interleukins,  heat  shock 
proteins  and  S 1 00  proteins.  None  of  the  pathways  showed  transient  regulation 
suggesting  that  the  tissue  response  is  sustained  past  the  24  time  point  in  this  data  set.  The 
response  of  the  gene  groups  over  time  indicate  that  tissue  repair  and  DNA  remodeling  are 
acute  and  may  be  completed  in  the  first  24  hours  post  exposure.  The  decrease  in  TNF 
and  interleukins  suggests  that  the  stress  response  is  dimished  while  the  tissue  attempts  to 
repair  DNA  and  other  damage.  A  response  of  this  type  is  considerably  different  that  what 
has  been  seen  in  high  dose  exposures,  and  supports  the  findings  of  other  recently 
emerging  in  vitro  data. 


Susanne  R.  Berglund,  Alison  Santana,  Dan  Li,  David  Rocke,  and  Zelanna  Goldberg, 
“Proteomic  Analysis  of  Low  Dose  Arsenic  and  Ionizing  Radiation  Exposure  on 
Keratinocytes,”  submitted  for  publication. 

The  human  health  effects  of  arsenic  and  ionizing  radiation  at  environmental  levels  have 
been  examined  separately,  but  little  information  regarding  potential  interactions  at  low 
doses  exist.  Arsenic  toxicity  may  be  synergistically  affected  by  concurrent  low  dose 
ionizing  radiation  especially  in  light  of  their  known  individual  carcinogenic  actions  at 
higher  doses.  A  proteomics  approach  was  used  to  examine  the  interaction  of  ionizing 
radiation  and  arsenic  in  a  human  cell  line.  Exposure  resulted  in  differential  proteomic 
expression.  Based  on  MSMS,  immunoblotting  and  statistical  analysis  (ANOVA); 
thirteen  proteins  were  identified  including:  annexin  XI,  calmodulin,  cyclophilin  A,  a- 
enolase,  epidermal  fatty  acid-binding  protein  (E-FABP),  heat  shock  protein  27,  histidine 
triad  nucleotide-binding  protein  1,  lactate  dehydrogenase  A  (LDH  A),  profilin-1,  protein 
disulfide  isomerase  precursor  (PDI),  pyruvate  kinase  M  isozyme,  R3372_l  and  S100A9. 
Six  of  the  proteins  responded  to  ionizing  radiation  and  eight  reacted  to  arsenic  exposure. 
Four  of  these  proteins,  a-enolase,  E-FABP,  LDH  A  and  PDI  had  an  effect  that  was 
different  than  would  be  expected  based  on  the  effects  of  either  toxicant  alone 
demonstrating  a  possible  response  to  the  combined  insult  that  is  substantially  different 
from  either  separate  treatment. 


Adaptive  Response  in  Epithelial  Tissue 


It  is  known  that  85%  of  cancers  originate  in  epithelial  tissues,  yet  there  is  little  literature 
on  adaptive  response  in  epithelial  cells.  The  most  common  cell  types  for  adaptive 
response  studies,  and  studies  of  low-dose  ionizing  radiation  in  general,  are  fibroblasts, 
lympocytes,  and  lymphoblastoid  cell  lines.  A  few  studies  have  been  conducted  on 
Chinese  hampster  ovary  (CHO)  cells,  which  are  fibroblast-like,  and  neuroblastoid  cell 
lines. 

Under  the  funding  from  this  project,  we  have  begun  investigations  into  adaptive  response 
in  epithelial  cells  (several  keratinocyte  cell  lines),  stromal  cells  (several  fibroblast  cell 
lines,  and  the  interactions  between  them.  We  have  used  the  MTS  assay,  a  micronuclei 
assay,  and  several  flow  cytometric  assays  to  refine  the  priming/challenge  dose  “sweet 
spot”  for  adaptive  response.  Although  the  data  are  not  yet  published,  we  have  good 
evidence  that  the  optimal  priming  dose  is  in  the  5-30  cGy  range  for  protection  against  a 
variety  of  challenge  doses,  depending  on  the  cell  type.  Illustrative  data  are  shown  in  the 
two  Figures  below. 


Adaptive  Response  with  4  Gy  Challenge  Dose 


Figure  1.  Micronuclei  assay  results  for  a  keratinocyte  cell  line  (HaCaTs)  at  priming  doses 
from  0  to  80  cGy  showing  that  the  optimal  priming  dose  is  probably  between  1  and  30 
cGy.  The  results  are  still  preliminary,  but  have  been  partially  confirmed  as  shown  in 
Figure  2. 
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Figure  2.  Micronuclei  assay  results  for  a  keratinocyte  cell  line  (HaCaTs)  at  priming  doses 
from  0  to  20  cGy  showing  that  the  optimal  priming  dose  is  probably  between  10  and  30 
cGy. 
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Human  In  vivo  Dose- Response  to  Controlled,  Low- Dose  Low  Linear 
Energy  Transfer  Ionizing  Radiation  Exposure 

Zelanna  Goldberg,1  David  M.  Rocke,2  Chad  Schwietert,1  Susanne  R.  Berglund,1  Alison  Santana,1 
Angela  Jones,1  Jorg  Lehmann,1  Robin  Stern,1  Ruixiao  Lu,3  and  Christine  Hartmann  Siantar4 


Abstract  Purpose:  The  effect  of  low  doses  of  low -linear  energy  transfer  (photon)  ionizing  radiation 
(LDIR,  <10  cGy)  on  human  tissue  when  exposure  is  under  normal  physiologic  conditions  is  of 
significant  interest  to  the  medical  and  scientific  community  in  therapeutic  and  other  contexts. 
Although,  to  date,  there  has  been  no  direct  assessment  of  the  response  of  human  tissue  to 
LDIR  when  exposure  is  under  normal  physiologic  conditions  of  intact  three-dimensional  archi¬ 
tecture,  vasculature,  and  cell-cell  contacts  (between  epithelial  cells  and  between  epithelial  and 
stromal  cells). 

Experimental  Design:  In  this  article,  we  present  the  first  data  on  the  response  of  human  tissue 
exposed  in  vivo  to  LDIR  with  precisely  controlled  and  calibrated  doses.  We  evaluated  transcrip- 
tomic  responses  to  a  single  exposure  of  LDIR  in  the  normal  skin  of  men  undergoing  therapeutic 
radiation  for  prostate  cancer  (research  protocol,  Health  Insurance  Portability  and  Accountability 
Act -compliant.  Institutional  Review  Board -approved).  Using  newly  developed  biostatistical 
tools  that  account  for  individual  splice  variants  and  the  expected  variability  of  temporal  response 
between  humans  even  when  the  outcome  is  measured  at  a  single  time,  we  show  a  dose-response 
pattern  in  gene  expression  in  a  number  of  pathways  and  gene  groups  that  are  biologically  plausi¬ 
ble  responses  to  LDIR. 

Results:  Examining  genes  and  pathways  identified  as  radiation-responsive  in  cell  culture  models, 
we  found  seven  gene  groups  and  five  pathways  that  were  altered  in  men  in  this  experiment. These 
included  the  Akt/phosphoinositide-3-kinase  pathway,  the  growth  factor  pathway,  the  stress/ 
apoptosis  pathway,  and  the  pathway  initiated  by  transforming  growth  factor- £  signaling,  whereas 
gene  groups  with  altered  expression  included  the  keratins,  the  zinc  finger  proteins  and  signaling 
molecules  in  the  mitogen-activated  protein  kinase  gene  group.  We  show  that  there  is  considerable 
individual  variability  in  radiation  response  that  makes  the  detection  of  effects  difficult,  but  still 
feasible  when  analyzed  according  to  gene  group  and  pathway. 

Conclusions:  These  results  show  for  the  first  time  that  low  doses  of  radiation  have  an  identifi¬ 
able  biosignature  in  human  tissue,  irradiated  in  vivo  with  normal  intact  three-dimensional  archi¬ 
tecture,  vascular  supply,  and  innervation.  The  genes  and  pathways  show  that  the  tissue  (a)  does 
detect  the  injury,  (b)  initiates  a  stress/inflammatory  response,  (c)  undergoes  DNA  remodeling, 
as  suggested  by  the  significant  increase  in  zinc  finger  protein  gene  expression,  and  (d)  initiates 
a  “pro-survival"  response.  The  ability  to  detect  a  distinct  radiation  response  pattern  following 
LDIR  exposure  has  important  implications  for  risk  assessment  in  both  therapeutic  and  national 
defense  contexts. 


Although  there  is  increasing  concern  and  interest  on  the  effects 
of  low-dose  low -linear  energy  transfer  ionizing  radiation 
(LDIR)  in  humans,  particularly  with  respect  to  secondary 
radiation  carcinogenesis  following  therapeutic  radiation,  there 
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exists  no  direct  evidence  that  doses  in  the  range  of  1  to  10  cGy 
have  any  biological  effects  whatsoever.  Studies  in  this  regard 
have  long  been  stymied  by  ethical  constraints  that  prohibit  the 
intentional  radiation  of  healthy  tissue  to  simply  examine  the 
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response  to  ionizing  radiation.  Attempts  to  address  this 
information  gap  have  approached  the  issue  from  two  different 
angles. 

Epidemiologic  studies  of  exposed  populations  have  evaluated 
global  population  exposure  and  health  outcomes  of  the  group, 
usually,  although  not  exclusively,  malignancy  (1-5).  A  careful 
review  of  these  studies  reveals  that,  at  the  lowest  exposure 
levels,  no  increase  in  detectable  risk  is  observed.  Whether  this  is 
a  true  threshold  or  simply  an  inherent  limitation  of  the 
epidemiologic  tools  is  the  subject  of  vigorous  debate. 

The  other  approach  to  determining  any  potential  biological 
activity  in  response  to  LD1R  has  been  through  classic  in  vitro  cell 
biology,  and  more  recently,  molecular  biology  using  genomics 
to  characterize  the  response  pattern.  Classic  end  points  such  as 
cell  lethality  have  not  been  sensitive  enough  to  show  any 
biological  effects  from  such  low  dose  exposures.  Using  genome¬ 
wide  techniques  on  cell  culture  models,  evidence  of  transcrip¬ 
tional  changes  following  LDIR  exposure  has  been  shown. 
However,  such  responses  are  often  subtle  and  inconsistent 
across  models  (6,  7). 

Cell  culture  data  of  primary  human  keratinocytes  and 
fibroblasts  exposed  to  low  or  high  doses  of  ionizing  radiation 
have  suggested  a  different  response  profile  under  each  of  these 
sets  of  circumstances  (6,  7).  Another  group  has  examined  brain 
samples  from  whole  animal  exposures  and  showed  a  tran- 
scriptomic  response  profile  to  LDIR  that  was  qualitatively 
different  than  that  seen  at  higher  doses  (8).  Each  of  these  studies 
has  identified  three  distinct  groups  of  transcripts:  those  that  were 
altered  by  high-  and  low-dose  exposures,  those  that  were  unique 
to  high-dose,  and  those  unique  to  the  low-dose  exposures.  This 
suggests  the  possibility  of  a  low-dose  radiation  response  that  is 
biologically  dissimilar  to  higher  doses  of  radiation. 

While  these  results  are  intriguing,  it  remains  unknown  how 
applicable  they  are  to  the  actual  human  response  because  the 
experiments  were  conducted  on  single  cell  layer,  single  cell  type 
in  vitro  models,  which  fail  to  accommodate  the  known 
importance  of  cell-cell  communication  in  radiation  response. 
Furthermore,  data  suggesting  that  transcriptional  changes  of 
key  response  genes,  e.g.,  protein  kinase  C,  were  markedly 
different  in  sensitivity  and  isoform  when  ionizing  radiation  was 
delivered  in  vitro  or  in  vivo  in  a  mouse  model.  This  highlights 
the  necessity  for  actual  human  data  to  begin  assessment  of  the 
human  responses  to  LDIR  (9). 

The  studies  described  herein  begin  to  address  this  informa¬ 
tion  gap.  We  have  developed  a  scientifically  valid  and  ethically 
compliant  model  for  direct  evaluation  of  LDIR  effects  in  normal 
human  tissue,  as  well  as  a  statistically  validated  tool  to  evaluate 
the  response.  Using  genomic  evaluation,  with  a  priori  selection 
of  transcripts  potentially  sensitive  to  radiation-induced  changes, 
we  have  (a)  established  the  first  data  set  of  human  responses  to 
well-defined  doses  of  LDIR,  ( b )  developed  a  new  methodology 
to  address  the  natural  heterogeneity  of  response  in  an  unselected 
patient  population,  (c)  developed  a  benchmark  data  set  of 
transcript  responses  to  LDIR  in  human  skin,  the  first  line  of 
defense  to  radiation  exposure  from  natural,  medical,  or  terrorist 
sources,  and  ( d )  showed  that  although  the  in  vitro  evaluation  of 
LDIR  effects  has  identified  genes  and  pathways  activated  by  such 
exposure,  it  does  not  characterize  the  response  in  tissue. 
Although  in  vitro  assays  do  serve  an  important  role  in  the 
support  of  translational  radiation  response  models,  they  cannot 
evaluate  true  radiation  response  in  humans. 


Materials  and  Methods 

Patients.  Men  from  the  clinical  practice  of  the  first  author,  who  are 
undergoing  therapeutic  radiation  for  the  treatment  of  their  early  stage 
prostate  cancer,  were  approached  for  possible  trial  participation. 
Informed  consent  was  obtained  (Institutional  Review  Board  -  approved, 
Health  Insurance  Portability  and  Accountability  Act -compliant).  At  the 
time  of  patient  treatment  planning  for  their  therapeutic  ionizing 
radiation,  an  extra  computed  tomography  scan  was  obtained  with  the 
patient  in  the  treatment  position.  This  scan  was  similar  to  the  one  used 
for  therapeutic  treatment  planning  purposes  but  was  without  any 
contrast  material  and  with  a  slab  of  tissue-equivalent  bolus  material 
over  the  areas  where  biopsies  were  to  be  taken.  The  research  scan  was 
used  for  protocol  treatment  planning.  The  bolus  material  on  the  skin 
ensured  that  the  region  was  in  electronic  equilibrium,  reducing  the 
radiation-absorbed  dose  uncertainty  in  the  sample.  The  biopsy  site  was 
determined  by  PEREGRINE  Monte  Carlo  dose  calculations  (10).  On  the 
first  day  of  treatment,  biopsy  points  were  identified,  and  marked  on 
the  patient  skin.  TLD  and/or  MOSFET  microdosimeters  were  placed  on 
the  marked  sites  to  confirm  the  accuracy  of  the  dose  delivered.  Biopsies 
were  done  and  the  samples  were  placed  in  RNAlater  (Ambion,  TX)  until 
processed.  Each  man  underwent  four  full  skin  thickness  punch  biopsies 
(3  mm  diameter  cores)  of  the  normal  abdominal  skin.  The  first  biopsy 
was  prior  to  any  ionizing  radiation  and  served  as  the  control  sample. 
The  other  three  samples  were  obtained  3  hours  following  ionizing 
radiation  at  sites  determined  to  have  received  1,  10,  or  100  cGy  point 
dose  at  the  skin  surface.  Biopsies  were  stored  at  -20°C  until  further 
processing. 

Tissue  disruption.  Each  biopsy  sample  was  loaded  into  a  Lysing 
Matrix  D  tube  (Qbiogene,  Irvine,  CA)  containing  1  mL  of  a  guanidine 
thiocyanate  solution  (5.1  mol/L  guanidine  thiocyanate,  50  mmol/L 
sodium  citrate,  50  mmol/L  EDTA,  and  0.5%  p-mercaptoethanol).  Each 
sample  underwent  three  rounds  of  bead  beating  in  the  Fastprep 
beadbeater  instrument  at  a  setting  of  6  m/s  (model  120A,  Qbiogene). 
Samples  were  kept  on  ice  between  each  pulverization  step.  After  tissue 
disruption,  the  samples  were  centrifuged  at  12,000  x  g  for  5  minutes  at 
4°C.  The  supernatant  was  removed  and  transferred  to  a  new  1.5  mL 
Eppendorf  tube. 

RNA  extraction  and  quantitation.  An  equal  volume  (1  mL)  of 
phenol  chloroform  (5:1,  pH  4.7)  containing  0.5%  n-lauroylsarcosine 
was  added  to  each  sample  and  incubated  for  5  minutes  at  room 
temperature  (Sigma,  St.  Louis,  MO).  The  samples  were  centrifuged  at 
12,000  x  g  for  10  minutes  at  2°C.  The  aqueous  phase  was  separated 
from  the  organic  phase  containing  DNA  and  protein.  An  additional 
phenol  chloroform  (5:1,  pH  4.7)  extraction  was  done  for  further 
purification.  One  milliliter  of  isopropanol  and  200  pL  of  3  mol/L 
sodium  acetate  was  added  to  each  sample  for  precipitation  and  stored 
overnight  at  -20°C.  The  precipitated  RNA  samples  were  cleaned  up 
using  RNeasy  columns  (Qiagen,  CA)  by  following  the  manufacturer's 
cleanup  procedure.  Modifications  to  the  procedure  included  passing  the 
RNA  through  the  binding  column  twice  to  increase  the  binding 
efficiency  and  elution  of  the  RNA  with  RNA  storage  solution  (Ambion). 
The  Turbo  DNA-Free  kit  was  used  to  remove  any  residual  DNA 
contamination  (Ambion).  RNA  integrity  was  verified  with 
the  Molecular  Devices  SpectraMax  plus.  The  RNA  samples  averaged 
1.6  pg  of  total  RNA  and  the  260:280  nm  ratios  were  between  1.7  and 
1.9.  The  RNA  was  stored  at  -80°C. 

Amplification  and  labeling.  The  RNA  was  processed  for  the 
Affymetrix  arrays  using  their  one-cycle  in  vitro  transcription  labeling 
kit  (Affymetrix,  Santa  Clara,  CA).  Approximately  500  pg  of  mRNA  was 
reverse-transcribed  using  a  T7-01igo(dT)  promoter  primer  in  the  first- 
strand  synthesis.  Second-strand  cDNA  synthesis  was  done  and  the 
resulting  double-stranded  cDNA  was  purified.  This  cDNA  served  as  a 
template  in  the  in  vitro  transcription  reaction.  The  in  vitro  transcription 
reaction  was  carried  out  in  the  presence  of  T7  RNA  polymerase  and  a 
biotinylated  nucleotide  analogue/ribonucleotide  mix.  The  biotinylated 
cRNA  targets  were  then  cleaned  up  and  fragmented  using  standard 
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Affymetrix  materials  and  protocols.  The  cRNA  was  hybridized  to  the 
Human  Genome  U133  Plus  2.0  arrays. 

Data  set  selection.  To  address  the  complexity  of  the  data  which 
would  be  generated  by  the  Affymetrix  GeneChips;  a  priori,  we  selected 
genes  and  pathways  to  be  examined.  These  genes  or  gene  pathways  were 
chosen  based  on  the  published  literature  in  radiation  response,  mostly, 
but  not  exclusively,  at  higher  doses.  To  address  the  variable  temporal 
response  between  individuals,  any  gene  in  which  there  was  evidence  in 
the  literature  to  justify  it  as  a  radiation-responsive  gene  was  evaluated, 
as  was  that  entire  gene  pathway.  This  was  to  capture  the  data  from  all  of 
the  sample  cohorts  who  were  evaluated  (biopsies  obtained)  at  a  fixed 
point  in  time  after  ionizing  radiation  exposure,  recognizing  that  the  rate 
of  the  transcriptional  response  would  vary  across  the  population. 
Below,  we  review  the  chosen  genes/pathways  and  their  rationales. 

Ionizing  radiation  has  been  shown  to  activate  nuclear  factor  kB 
through  the  degradation  of  inhibitor-KB.  Therefore,  the  signaling 
pathways  upstream  of  this  (extracellular  signal-regulated  kinase  and 
mitogen-activated  protein  kinase)  were  also  considered  in  the  pathways 
analysis  (11).  Cell  cycle- modulating  genes  such  as  cyclins,  cdc25 
phosphatases,  GADD45,  ATM,  and  ATR  as  well  as  the  chk-1  and  chk-2 
transcripts  were  also  examined.  The  TP53  pathway,  with  upstream  and 
downstream  signaling  and  effector  transcripts,  was  examined,  as  was 
transforming  growth  factor-p  and  the  ubiquitons  (12,  13).  Growth 
factors,  such  as  basic  fibroblast  growth  factor,  and  its  downstream 
targets  (RAS/MEK/mitogen-activated  protein  kinase/RSK)  were  assessed, 
as  well  as  other  apoptosis-related  paths  and  survival  signaling  pathways 
(Akt/phosphoinositide-3-kinase,  growth  factors  such  as  tumor  necrosis 
factor,  epidermal  growth  factor  receptor,  vascular  endothelial  growth 
factor,  and  their  downstream  effectors  Ras,  Raf-ERK)  were  evaluated 
(14-18).  DNA  repair  genes,  in  particular,  the  Rad  52  epistasis  group, 
have  been  reported  by  some  to  be  radiation-responsive,  although  others 
have  questioned  this  relationship  given  the  constituent  expression  of 
DNA  repair  proteins  (6,  19).  Inflammatory  mediators  such  as 
cyclooxygenase-2,  prostaglandin  E2,  and  interleukins  have  also  been 
reported  (20).  We  also  selected  some  tissue-specific  genes  with 
structural  functions  (i.e.,  keratins)  given  that  we  anticipated  some 
tissue  reorganization  in  response  to  ionizing  radiation.  This  type  of 
response  was  reported  by  Ding  et  al.  with  increased  expression  of 
cytoskeleton  components,  ANLN  and  KRT15,  from  normal  human 
fibroblasts  (6).  DNA  damage/remodeling  via  topoisomerase  and  zinc 
finger  proteins  has  also  been  reported  and  hence  these  transcripts  were 
also  examined  (21)  along  with  genes  involved  in  cell-cell  signaling, 
such  as  GRAP2  and  GPR51,  and  cytoskeleton,  such  as  ANLN  and 
KRT15. 

Therefore,  there  were  seven  preidentified  gene  pathways  analyzed. 
They  include  the  Akt/phosphoinositide-3-kinase  pathway  (14,  18,  22,  23), 
chemokine  pathway  (9,  24),  fibronectin  pathway  (14,  24,  25),  growth 
factor/insulin  pathway  (9),  inflammation  pathway  (9,  14,  20, 
22,  24,  26),  stress/apoptosis  pathway  (9,  13,  14,  20,  22-24,  27),  and 
the  transforming  growth  facto r-p/cyclin/ubiquitin  pathway  (11,  13,  14, 
18,  22,  24,  25,  28).  There  were  19  gene  groups  analyzed  (see  Table  1). 

Statistical  methods.  The  data  for  this  analysis  consist  of  32 
Affymetrix  HGU133  Plus  2.0  GeneChips,  four  from  each  of  eight 
patients,  at  doses  of  1,  10,  and  100  cGy,  as  well  as  a  preexposure  control 
at  0  cGy.  The  HGU133  Plus  2.0  contains  1,354,896  probes  divided  into 
54,675  probe  sets,  which  have  been  summarized  using  the  GLA 
expression  index  (29),  although  the  results  are  very  similar  if  RMA  is 
used  instead  (30). 

The  analysis  looked  for  dose-response  patterns  that  were  linear  in 
dose  either  on  the  original  cGy  scale,  or  using  a  modified  log  dose 
(MLD)  in  which  the  positive  exposures  are  coded  as  0,  1,  and  2,  and  the 
0  exposure  is  coded  as  -1,  as  if  it  were  0.1  cGy.  Because  we  expected 
(and  found)  variability  between  individuals  as  to  timing  of  response  on 
pathways  and  on  the  type  of  transcript  within  classes  of  genes,  we  did 
an  analysis  in  which  we  aggregated  results  across  individuals  and 
transcripts  in  order  to  amplify  possibly  weak  signals,  without  inducing 
false-positives. 


Full  details  are  given  in  Rocke  (31)  and  Rocke  et  al.  (32),  but  in 
summary,  for  each  probe  set  and  each  individual,  we  conducted  a  linear 
regression  and  computed  the  t  statistic  testing  for  whether  there  is  up- 
regulation,  down-regulation,  or  neither.  The  t  score  will  be  positive  if 
there  is  a  trend  towards  up-regulation,  and  negative  if  there  is  an 
opposite  trend.  We  examined  the  entire  collection  of  t  scores  for  a  given 
gene  group  or  pathway  to  determine  the  possibility  of  a  weak  signal  in 
the  direction  of  up-regulation  or  in  the  direction  of  down-regulation. 
We  did  this  by  testing  the  hypothesis  that  the  collection  of  t  scores  had 
a  mean  of  zero  using  the  one-sample  t  test  or  a  median  of  zero  using  the 
Wilcoxon  rank-sum  test.  If  this  was  rejected,  it  indicated  that  there  was  a 
trend  towards  up-regulation  with  dose  of  the  gene  group  or  pathway  if 
the  t  scores  are  biased  in  a  positive  direction  or  towards  down- 
regulation  if  the  t  scores  are  biased  in  a  negative  direction.  Because  of 
possible  correlations  in  these  test  statistics,  the  most  reliable  P  values 
are  provided  by  resampling  by  repeatedly  sampling  random  groups  of 
transcripts  of  the  same  size.  The  empirical  P  value  is  the  fraction  of  cases 
in  the  resampled  gene  sets  in  which  the  test  statistic  was  more 
significant  than  the  actual  test  statistic  from  the  gene  group  or  pathway. 
These  were  always  two-sided  P  values. 

Results 

The  Affymetrix  HGU133  Plus  2.0  array  platform  shows  the 
advantages  and  disadvantages  of  whole  genome  assays.  The 
greatest  advantage  of  genome-wide  assays  is  the  comprehen¬ 
siveness.  Not  only  are  the  genes  represented  on  the  Affymetrix 
chip,  but  each  gene  is  represented  by  multiple  exons.  The 
disadvantage  is  that  with  such  a  large  number  of  probe  sets, 
there  will  be  many  apparently  significant  changes  that  occur 
strictly  by  chance,  and  to  control  for  this  requires  a  very  high 
bar  for  declaring  differential  expression  to  be  significant.  The 
model  that  uses  all  32  expression  values  for  each  probe  set  and 
tests  for  significant  patient  effects  and  dose-response  effects 
would  perhaps  be  sensitive  enough  if  the  response  were  highly 
consistent  in  time  across  individuals,  but  if  the  response  is 
varied  or  diffuse,  this  may  not  yield  significant  results.  In  the 
54,765  probe  sets,  and  using  the  methods  of  Rocke  (31),  we 
found  no  cases  in  which  the  overall  dose-response  coefficient  is 
significant  at  a  5%  FDR  level,  and  1 1  where  the  FDR-adjusted 
significance  level  is  better  than  10%  (31,  33,  34).  The  patient 
effect,  however,  was  significant  at  the  5%  FDR-adjusted  level  for 
13,514  of  the  54,675  probe  sets,  and  significant  at  the  10% 
FDR-adjusted  level  for  18,605  of  the  probe  sets.  This  shows 
how  important  the  individual  variation  is  in  gene  expression  in 
radiation  response.  Thus,  we  chose  to  analyze  the  results  using 
the  methods  of  Rocke  et  al.  (32)  and  using  the  gene  groups  and 
pathways  identified  in  the  literature,  as  described  in  Materials 
and  Methods  above  (32). 

Tables  1  and  2  show  the  results  for  the  preidentified  gene 
groups  and  pathways.  Of  the  19  gene  groups,  7  have 
empirical  P  <  0.0125  on  either  the  regression  on  dose  or 
the  regression  on  MLD,  using  either  the  t  test  or  the  Wilcoxon 
test,  a  number  far  exceeding  chance  (the  expected  number  out 
of  19  is  <1,  and  the  probability  of  getting  7  significant  by 
chance  is  <0.0003).  The  largest  effects  are  for  the  BCL-6 
group,  keratins,  protein  disulfide  isomerases,  SI 00,  and  zinc 
finger  proteins. 

For  the  pathways,  five  of  seven  pathways  show  significant 
difference  at  the  0.0125  level  in  dose  or  MLD  using  either  the 
t  test  or  Wilcoxon  test,  with  an  expected  number  <1,  and  the 
probability  of  having  five  or  more  significant  being 
-0.00006.  The  effect  is  strongly  exemplified  by  transforming 
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Table  1.  Results  for  preidentified  gene  groups 

Gene  group 

No.  of 

Direction  of 

ToTS  t.  P  value 

ToTS  W, 

Empirical  t. 

Empirical  W, 

probe  sets 

effect 

(dose/MLD) 

P  (dose/MLD) 

P  (dose/MLD) 

P  (dose/MLD) 

BAG 

9 

0.1480 

0.0397 

0.1865 

0.0460 

0.0251 

0.0332 

0.0195 

0.0420 

BCL2 

47 

0.2389 

0.6906 

0.3410 

0.7400 

0.6898 

0.8323 

0.7160 

0.9960 

BCL6 

8 

Up 

0.0005 

0.0005 

0.0005 

0.0000 

0.0091 

0.0039 

0.0035 

0.0020 

BMPs 

7 

0.2181 

0.2227 

0.2705 

0.2360 

0.5623 

0.5434 

0.6275 

0.5440 

BMP  receptor 

8 

0.4613 

0.7660 

0.5320 

0.7360 

0.6515 

0.3828 

0.6815 

0.4480 

Cyclins 

123 

0.1200 

0.3530 

0.2370 

0.3400 

0.1347 

0.2394 

0.1165 

0.3560 

Cytokines 

86 

Up 

0.0169 

0.0013 

0.0365 

0.0020 

0.0184 

0.0074 

0.0140 

0.0260 

GADD  45 

6 

0.0563 

0.0668 

0.0555 

0.0820 

0.0247 

0.0380 

0.0175 

0.0500 

Heat  shock 

58 

0.7280 

0.6620 

0.7725 

0.7100 

0.5784 

0.3377 

0.6030 

0.2740 

Interleukin 

147 

0.9877 

0.4346 

0.9935 

0.5540 

0.6330 

0.9583 

0.6395 

0.8080 

Keratins 

101 

Down 

0.0020 

7.21  E-1 5 

0.0070 

0.0000 

2.99E-14 

2.20E-16 

0.0000 

0.0000 

Mitogen-activated  protein  kinase 

131 

Up 

0.0020 

0.0025 

0.0105 

0.0040 

0.0143 

0.0164 

0.0140 

0.0440 

Protein  disulfide  isomerase 

8 

Down 

0.0139 

0.0240 

0.0150 

0.0420 

0.0121 

0.0064 

0.0070 

0.0000 

RAD  51 

10 

0.4673 

0.2794 

0.5340 

0.3220 

0.0411 

0.1495 

0.0365 

0.1900 

S100 

21 

Down 

0.0144 

0.0015 

0.0020 

0.0002 

0.0001 

0.0001 

0.0000 

0.0000 

Serine/threonine  kinase 

73 

0.3329 

0.9555 

0.4425 

0.9340 

0.3778 

0.9997 

0.4090 

0.8860 

Tumor  necrosis  factor 

109 

0.8158 

0.2513 

0.8710 

0.3520 

0.8300 

0.8376 

0.8505 

0.6220 

Topoisomerase 

14 

0.2022 

0.1399 

0.2490 

0.1960 

0.0541 

0.1561 

0.0500 

0.2520 

Zinc  finger 

799 

Up 

2.68E-05 

7.79  E- 5 

0.0135 

0.0000 

4.17E-06 

2.20E-16 

0.0000 

0.0000 

NOTE:  The  empirical  P  values  are  based  on  simulations  of  2,000  trials  for  the  t  test  and  1,000  trials  for  theWilcoxon  (W)  test.  Rows  in  bold  are  gene  groups  for  which 
there  is  good  evidence  of  differential  expression.  Other  gene  groups  with  borderline  significance  include  BAG,  GADD  45,  RAD51,  and  topoisomerase.  In  each  cell,  the 

upper  P  value  is  for  regression  on  dose,  and  the  lower  P  value  for  the  regression  on  MLD. 

growth  factor-p/cyclin/ubiquitin  pathway  activated  in  re¬ 
sponse  to  stress/apoptosis,  inflammation,  growth  factor/ 
insulin,  and  Akt/phosphoinositide-3-kinase  pathways.  Of  the 
seven  preidentified  gene  pathways  analyzed,  five  showed  a 
significant  dose-response  up-regulation  to  LDIR  exposure.  The 
two  that  did  not  were  the  chemokine  and  fibronectin 
pathways. 

There  were  1 9  gene  groups  analyzed,  of  which  seven  showed 
radiation-responsive  differential  gene  expression.  The  groups 
that  were  up-regulated  include  BCL6  (14),  cytokines  (9,  23,  25, 
28),  mitogen-activated  protein  kinases  (14,  20,  23,  24,  35),  and 


zinc  finger  proteins  (27).  Those  that  were  down-regulated 
include  the  keratins  (25),  the  protein  disulfide  isomerases  (15), 
and  S100  (36).  Those  that  did  not  change  significantly 
included  BAG  (37),  BCL2  (38),  BMPs  (39,  40),  BMP  receptors 
(40),  cydins  (41,  42),  GADD45  (43),  heat  shock  (41,  44), 
interleukins  (45),  RAD51  (19),  serine/threonine  kinases 
(46,  47),  tumor  necrosis  factor-p  (48),  and  topoisomerases 
(21,  49).  The  radiation-responsive  gene  groups  and  pathways 
identified  in  our  human  biopsy  studies  have  key  regulatory 
roles  in  cell  cycle  arrest,  DNA  repair,  inflammation,  and 
apoptosis. 
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Table  2.  Results  for  preidentified  pathways 

Pathway 

No.  of 

Direction  of 

ToTS  t 

ToTS  W, 

Empirical  t. 

Empirical  W, 

probe  sets 

effect 

P  (dose/MLD) 

P  (dose/MLD) 

P  (dose/MLD) 

P  (dose/MLD) 

Akt/phosphoinositide-3-kinase  pathway 

99 

Up 

0.0014 

0.0002 

0.0050 

0.0020 

0.1333 

0.0006 

0.1335 

0.0020 

Chemokine  pathway 

79 

0.2716 

0.4047 

0.3895 

0.4780 

0.8565 

0.2649 

0.8730 

0.4080 

Fibronectin  pathway 

196 

0.0206 

0.0649 

0.0990 

0.1020 

0.0497 

0.0169 

0.0455 

0.0460 

Growth  factor/insulin  pathway 

208 

Up 

0.0001 

0.0072 

0.0015 

0.0100 

0.5044 

0.0083 

0.5405 

0.0220 

Inflammation  pathway 

78 

Up 

0.0365 

0.0013 

0.0785 

0.0020 

0.0687 

0.1187 

0.0810 

0.1780 

Stress/apoptosis  pathway 

151 

Up 

0.9343 

0.0739 

0.9500 

0.0400 

0.0359 

0.0029 

0.0420 

0.0070 

Transforming  growth  factor-fi/ 

355 

Up 

0.0003 

0.0005 

0.0095 

0.0000 

cyclin/ubiquitin  pathway 

0.0014 

0.0012 

0.0005 

0.0010 

NOTE  The  empirical  P  values  are  based  on  simulations  of  2,000  trials  for  the  t  test  and  1,000  trials  for  the  Wilcoxon  (W)  test.  Rows  in  bold  are  pathways  for  which  there  is 

good  evidence  of  differential  expression.  In  each  cell,  the  upper  P  value  is  for  regression  on  dose,  and  the  lower  P  value  for  the  regression  on  MLD. 

In  most  cases,  gene  groups  and  pathways  that  were  significant 
over  all  patients  were  also  significant  for  individual  patients, 
although  this  did  not  always  occur,  which  is  consistent  with  the 
assumption  of  interindividual  variability  in  response.  In  many 
cases,  the  individual  tests  were  not  significant,  consistent  with 
the  lower  statistical  power  of  a  test  based  on  one  individual,  and 
in  only  a  few  cases  was  the  effect  significant  in  the  opposite 
direction.  This  overall  consistency  gives  additional  confidence 
in  the  realism  of  the  results.  Details  are  given  in  Table  3. 

As  noted  above  in  Materials  and  Methods,  all  expression  data 
were  analyzed  using  a  model  of  linear  dose-response.  Thus, 
significant  results  indicate  linear  response  in  expression  up  (or 
down)  with  dose  over  the  range  of  1  to  100  cGy.  Those  gene 
groups  or  pathways  that  were  not  statistically  significant  on  a 


linear  model  had  variable  responses  across  dose  and  could  not 
be  modeled  by  a  linear  fit. 

The  microarray  data  from  this  study  are  available  in 
Minimum  Information  About  a  Microarray  Gene  Experi¬ 
ment-compliant  format  from  the  National  Center  for 
Biotechnology  Information  Gene  Expression  Omnibus  repos¬ 
itory  and  include  the  32.  CEL  files  and  the  probe  set  summary 
data  for  all  32  arrays  and  54,675  probe  sets  in  a  Microsoft 
Excel  file,  as  well  as  the  experimental  metadata.  The  data  sets 
defining  the  gene  groups  and  pathways,  and  the  programs  in 
the  R  language  that  were  used  to  process  the  data  are  available 
as  supplementary  files  from  the  Journal  web  site.  All  of  these 
files  are  also  available  on  the  second  authors  web  site,  http:// 
www. idav.ucdavis.edu/~dmrocke. 


Table  3.  Results  for  individual  patients  for  significant  gene  groups  and  pathways 


Gene  group  or  pathway 

Direction  of 

effect 

No.  of  patients 
up-regulated 

No.  of  patients 
down-regulated 

No.  of  patients 
not  significant 

BCL6 

Up 

4 

0 

4 

Cytokines 

Up 

3 

1 

4 

Keratins 

Down 

8 

0 

3 

Mitogen-activated  protein  kinase 

Up 

2 

0 

6 

Protein  disulfide  isomerase 

Down 

0 

2 

6 

SI  00 

Down 

0 

3 

5 

Zinc  finger 

Up 

6 

1 

1 

Akt/phosphoinositide-3 -kinase  pathway 

Up 

5 

0 

3 

Growth  factor/insulin  pathway 

Up 

4 

0 

4 

Inflammation  pathway 

Up 

4 

1 

3 

Stress/apoptosis  pathway 

Up 

3 

1 

4 

Transforming  growth 

Up 

3 

0 

5 

factor-  p/cyclin/ubiquitin  pathway 
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Discussion 

Ionizing  radiation  is  a  naturally  ubiquitous  toxicant.  Back¬ 
ground  radiation  is  now  supplemented  through  most  peoples' 
lives  by  industrial  exposures  and  medical  exposures  for 
diagnostic  and/or  therapeutic  reasons.  Furthermore,  in  the  era 
of  intensity-modulated  radiation  therapy,  larger  volumes  of 
normal  tissue  are  exposed  to  low  doses  of  ionizing  radiation. 
There  is  also  a  newly  elevated  threat  from  radiologic  terrorism, 
driving  the  need  to  determine  the  biological  effects,  if  any,  of 
low  levels  of  radiation  exposure  in  the  human. 

By  devising  a  dosimetrically  sound  sampling  method  on 
individuals  who  are  otherwise  healthy  but  receiving  therapeutic 
localized  radiation  for  early  stage  prostate  cancer,  we  have 
developed  a  new  model  system  whereby  direct  evaluation  of 
radiation  effects  can  be  determined.  This  model  is  ethically 
sound,  Institutional  Review  Board  -  approved,  and  allows  real¬ 
time  sampling  of  human  tissue  after  in  vivo  radiation  exposures 
of  precisely  calibrated  radiation  doses.  This  model  is  therefore 
unique  in  allowing  the  evaluation  of  whole  tissue  effects  when 
exposure  is  under  normal  physiologic  conditions. 

Biopsies  were  obtained  from  patients  3  hours  after  an  acute 
exposure  to  LDIR.  This  time  point  was  chosen  for  practical 
considerations.  More  recently,  emerging  data  have  shown  that 
low  dose-response  may  be  better  defined  at  a  later  time  point 
(6,  7,  35).  Although  none  of  the  published  data  have  examined 
an  in  vivo  human  model,  it  is  reasonable  to  question  whether  a 
greater  response  may  be  seen  after  a  longer  interval.  If  so,  this 
would  suggest  that  secondary  tissue  effects,  i.e.,  cellular 
responses  to  signaling  from  neighboring  cells  may  be  as 
important  in  tissue  response  to  LDIR  as  the  initial  radiation 
deposition  within  the  cells.  Studies  to  better  characterize  the 
temporal  response  to  LDIR  in  tissue,  in  a  separate  cohort  of 
men,  are  ongoing  and  will  be  reported  separately. 

Using  genome-wide  analysis,  preselected  targeted  genes,  and 
pathways  from  in  vitro  studies,  we  have  shown  that  human  skin 
initiates  a  transcript  program  to  enhance  survival  whereas 
simultaneously  undergoing  an  inflammatory  and  stress  re¬ 
sponse.  Whether  these  pathways  are  all  activated  in  the  same 
cells  within  the  tissue,  or  represent  responses  from  different 
tissues  (i.e.,  epithelial  versus  stromal  versus  vasculature)  cannot 
be  determined  from  this  data  set.  Laser  capture  dissection  of  the 
biopsied  tissue  for  separate  analysis  of  each  cell  type  can  be 
undertaken  now  that  benchmark  data  have  been  compiled 
showing  that  low-dose  ionizing  radiation  is  biologically  active 
in  a  definable  way  in  the  human.  These  data  were  analyzed  on 
the  assumption  of  a  linear  dose-response.  Nonlinearity  at  the 
lowest  doses  cannot  be  excluded. 

The  tissues  sampled  and  examined  were  evaluated  as  whole 
tissue.  Therefore,  the  genomic  signature  is  a  summative 
evaluation  of  the  responses  of  thousands  of  cells  of  several 
lineages,  containing  both  epithelial  and  stromal  cells.  Some  of 
the  cells  would  represent  lineages  that  are  classically  "radiation- 
sensitive"  whereas  others  were  radiation-resistant.  Although 
this  model  does  not  allow  precise  mechanistic  intervention,  it 
does  offer  a  comparative  advantage  in  giving  us  tissue  level 
response  which  will  be  the  defining  level  of  human  response — 
i.e.,  the  tissue  level  response  is  more  than  just  the  summative 
whole  of  individual  cell  responses.  The  doses  examined  in  this 
study,  and  those  of  interest  to  the  radiation  risk  community  are 
well  below  those  that  are  frankly  cytotoxic.  Therefore,  the  tissue 


response  will  inevitably  be  an  amalgamation  of  the  primary 
response  of  the  cell  and  the  modulating  effects  of  the  cells 
surrounding  it,  both  of  the  same  lineage  and  others.  This 
amalgamated  response  therefore  incorporates  the  bystander 
response  and  represents  the  organismal  response  to  the 
radiation  insult.  Aside  from  better  determining  the  extent  of 
true  human  response  to  low  doses  of  ionizing  radiation,  it  also 
provides  a  potential  biosignature  for  evaluating  skin  exposure 
to  low  doses  of  ionizing  radiation,  well  below  doses  expected  to 
produce  symptomatic  radiation  effects. 

We  have  devised  a  methodology  whereby  an  individual 
biopsy  from  these  same  sample  types  can  be  processed  for  both 
RNA  and  protein  with  both  collected  for  analysis.  This  will 
allow  a  matched  genomic  and  proteomic  evaluation. s  These 
evaluations  are  ongoing,  but  will  further  extend  our  ability  to 
define  the  tissue  level  response  of  human  skin  to  low  doses  of 
ionizing  radiation. 

Defining  the  effects  of  low-dose  radiation  in  humans  requires 
a  model  with  diverse  genetic  background  with  the  subsequent 
statistical  tool  development  to  identify  the  effect  of  the  agent 
above  the  background  diversity  of  the  study  population. 
Therefore,  we  did  these  benchmark  studies  on  an  otherwise 
unselected  population  and  devised  statistically  rigorous  tools  to 
address  human  complexity.  Limitations  of  a  sample  set  consist¬ 
ing  of  males  with  early  stage  prostate  cancer  are  a  clear  lack  of 
gender  diversity  and  a  limited  age  range.  Radiation  sensitivities 
in  other  populations,  specifically  the  young  and  fetus'  in  utero , 
cannot  be  extrapolated  from  these  data.  This  data  set  included 
samples  from  six  Caucasians  and  two  African-Americans. 

The  skin  forms  the  first  barrier  of  the  body  to  radiation 
exposure  and  is  an  easy  tissue  to  biopsy  for  scientific 
evaluation,  clinical  studies  and  potential  field  application  for 
population  screening  for  radiation  exposure.  It  is  a  dynamic, 
complex  system  with  multiple  cell  types.  Because  only  healthy 
skin  was  biopsied  (not  previous  scars  or  cancerous  lesions),  this 
radiation  profile  may  represent  a  generalizable  biosignature. 
The  age  distribution  within  the  patient  population  sampled 
does  not  allow  an  analysis  of  age  effects,  and,  as  noted  above, 
gender  effects  cannot  be  excluded.  Relating  this  transcriptomic 
signature  to  a  functional  tissue  end  point  for  human  health  is 
part  of  an  ongoing  project. 

These  benchmark  data  show  for  the  first  time  that  low  doses 
of  radiation  have  an  identifiable  biosignature  in  human  tissue, 
irradiated  in  vivo  with  normal  intact  three-dimensional 
architecture,  vascular  supply,  and  innervation.  The  genes  and 
pathways  show  that  the  tissue  (a)  deteas  the  injury,  ( b ) 
initiates  a  stress/inflammatory  response,  (c)  undergoes  DNA 
remodeling,  suggested  by  the  significant  increase  in  genes  for 
zinc  finger  proteins  being  expressed,  and  ( d )  initiates  a  "pro- 
survival"  response. 
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Abbreviations: 


As 

sodium  arsenite 

CaM 

calmodulin 

E-FABP 

epidermal  fatty  acid-binding  protein 

HINT-1 

histidine  triad  nucleotide-binding  protein  1 

HSP27 

heat  shock  protein  27 

IR 

ionizing  radiation 

LDH-A 

lactate  dehydrogenase  A 

PDI 

protein  disulfide  isomerase  precursor 
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Summary 


While  the  human  health  effects  of  arsenic  and  ionizing  radiation  have  been  examined 
separately,  there  is  little  information  regarding  their  interactions  at  low  doses  approaching 
environmental  levels.  The  idea  of  synergistic  toxicity  is  a  well-known  phenomenon  in  human 
disease.  Arsenic  toxicity  may  be  affected  at  low  dose  exposure  by  concurrent  ionizing  radiation 
especially  in  light  of  their  known  individual  carcinogenic  actions  at  higher  doses. 

We  employed  a  proteomics  approach  in  a  human  keratinocyte  cell  line  to  examine  the 
interaction  of  ionizing  radiation  and  arsenic.  Toxicant  exposure  resulted  in  differential  proteomic 
expression  based  on  immunoblotting  and  statistical  analysis  (ANOVA)  of  identified  proteins. 

Thirteen  proteins  were  identified  including:  annexin  XI,  calmodulin,  cyclophilin  A,  a- 
enolase,  epidermal  fatty  acid-binding  protein,  heat  shock  protein  27,  histidine  triad  nucleotide¬ 
binding  protein  1,  lactate  dehydrogenase  A,  profilin-1,  protein  disulfide  isomerase  precursor, 
pyruvate  kinase  M  isozyme,  R3372_l  and  S100A9. 

Six  proteins  responded  to  ionizing  radiation  and  eight  reacted  to  arsenic  exposure.  Four 
proteins,  a-enolase,  epidermal  fatty  acid  binding  protein,  lactate  dehydrogenase  A  and  protein 
disulfide  isomerase,  had  an  interaction  effect  that  was  different  than  would  be  expected  based  on 
the  effects  of  either  toxicant  alone.  These  data  demonstrate  a  possible  response  to  the  combined 
insult  that  is  substantially  different  from  either  separate  treatment.  The  response  of  a-enolase 
and  protein  disulfide  isomerase  were  dissimilar  from  what  has  been  found  using  higher  doses, 
adding  to  the  idea  that  cellular  response  to  low  dose  exposures  are  distinct  from  the  responses  to 
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high  dose  exposures.  This  study  begins  to  address  the  cellular  proteomic  response  to  low  dose 
exposures  of  dual  environmental  toxicants. 

Introduction 

Low  level  arsenic  and  low  dose  ionizing  radiation  are  both  environmental  toxicants. 

While  data  exist  which  examine  the  human  health  effects  of  either  toxicant  separately,  there  are 
no  data  in  the  literature  regarding  possible  interactions  at  low  doses.  Yet,  interactions  of 
toxicants  to  produce  synergistic  toxicity  is  a  well-known  phenomenon  in  human  disease  (e.g. 
radon/tobacco  or  arsenic/tobacco  smoke  in  lung  cancer  induction)  ’  .  Additionally,  the  FDA  has 
approved  the  use  of  low  levels  of  arsenic  to  treat  acute  promyelocytic  leukemia.  With  the 
growing  use  of  IMRT,  which  results  in  an  increase  in  healthy  tissue  exposed  to  low  dose  ionizing 
radiation,  the  possibility  of  receiving  a  combined  or  sequential  exposure  to  these  particular 
toxicants  in  a  medical  setting  requires  further  biologic  characterization 3.  Arsenic  toxicity  may 
be  affected  at  low  dose  exposure  levels  by  concurrent  ionizing  radiation  especially  in  light  of 
their  known  carcinogenic  actions  individually  at  higher  doses.  Prior  studies  that  examine  each 
substance  independently  suggest  the  possibility  of  overlapping  cellular  defense  pathways  which 
are  activated  by  the  toxicants  4'7. 

Arsenic  is  a  naturally  occurring  metalloid  that  can  be  solubilized  in  water  under  certain 

O  Q 

conditions,  posing  the  threat  of  contaminated  drinking  water  ’  .  It  has  been  well  documented 

O  1  A  1  O 

that  exposure  to  arsenic  can  contribute  to  skin,  bladder,  liver  and  lung  cancers  ’  '  .  While  the 
mechanisms  of  arsenic  toxicity  are  still  not  fully  understood,  several  ideas  have  been  postulated 
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including  induction  of  oxidative  stress,  decreased  functioning  of  DNA  repair  systems, 
chromosomal  abnormality  and  altered  growth  factors  14'16.  There  is  evidence  that  exposure  to 
arsenic  generates  an  oxidative  stress  response  which  can  cause  an  increase  in  the  presence  of 
reactive  oxygen  species  (ROS) l3, 14.  Exposure  of  cells  to  ROS  can  have  many  effects,  including 
DNA  damage  and  mutation  17.  Changes  in  DNA  repair  may  occur  by  changing  the  expression  of 
genes  involved  in  the  synthesis  of  DNA  repair  enzymes,  rendering  cells  unable  to  fix  damaging 
mutations  ’  .  This  mechanism  has  also  been  proposed  as  an  explanation  to  arsenic’s  co- 
mutagenic  effects  l9. 

Ionizing  radiation  (IR)  exposure  is  unavoidable  in  the  environment  and  further  exposure 
is  commonly  obtained  through  medical  imaging.  Average  background  radiation  in  the  United 
States  is  360  mrem  (whole  body  equivalent  exposure)  with  medical  imaging  exposures  adding  up 
to  5  cGy  focal  absorbed  dose.  There  is  substantial  debate  regarding  the  biological  effects  of  low 
dose  IR  in  humans  and  no  direct  information  on  how  such  exposure  may  alter  the  response  of 
cells  to  other  environmental  toxicants  20'22.  As  with  arsenic,  IR  exerts  the  majority  of  its  toxicity 
through  the  intracellular  generation  of  ROS.  The  ROS  affect  the  target  cell,  as  well  as  generating 
a  redox  sensitive  signal  to  surrounding  cells,  the  so-called  bystander  effect  .  It  is  therefore 
plausible  that  oxidative  stress  induced  from  IR  could  substantially  enhance  the  effects  of 
otherwise  minimally  toxic,  sublethal  exposures  of  arsenic. 

Much  remains  uncertain  about  the  effects  of  these  toxicants  at  low  doses.  Historically, 
most  radiation  studies  have  involved  high  dose  exposures  with  the  assumption  that  the  toxicant 
profile  could  be  extrapolated  down  in  a  linear  manner  for  low  dose  exposures.  This  is  the 
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underlying  assumption  in  the  linear  no-threshold  model  of  radiation  effects,  which  is  currently 
undergoing  challenge  (as  reviewed  in: 24'27).  Doubt  still  remains  regarding  the  biological  effects 

of  low  dose  IR  on  humans,  how  the  effects  are  exerted  and  the  shape  of  the  low  dose  IR  response 

22 

curve 

As  with  IR  exposure,  much  remains  unknown  about  the  dose  response  curve  for  low  level 
arsenic  exposure.  In  the  past,  it  was  assumed  that  the  response  for  high  doses  of  arsenic  could 
also  be  applied  to  low  dose  exposures  9.  However,  current  studies  have  yet  to  demonstrate  a 
direct  relationship  between  low  dose  arsenic  exposure  and  cancer,  suggesting  a  nonlinear 
relationship  and  supporting  the  idea  that  the  dose  response  curve  at  low  levels  can  not  be  inferred 
from  high  dose  studies  l3.  Thus,  as  there  is  no  accepted,  comprehensive  model  describing  the 
mechanism  by  which  low  dose  toxin  exposures  exert  their  effects,  there  is  no  predictive 
modeling  that  can  address  the  potential  interacting  effects  of  co-exposures  on  cells.  Therefore, 
direct  empiric  study  remains  the  cornerstone  of  understanding  potential  interactions  and  how  the 
co-exposures  may  alter  the  safety  profile  of  each. 

The  skin  is  the  major  barrier  to  many  environmental  toxicants,  with  keratinocytes  being 
the  most  prevalent  cell  type.  Arsenic  exposure  can  lead  to  many  diseases  of  the  skin,  while  IR 
exposures  generally  penetrate  the  skin.  Chronic  arsenic  exposure  often  leads  to 
hyperpigmentation,  hyperkeratosis  and  arsenic  induced  Bowen’s  disease  15.  In  time,  Bowen’s 
disease  can  progress  into  invasive  skin  cancer  in  the  form  of  basal  or  squamous  cell  carcinoma 
15.  Keratinocytes  are  readily  cultivated  and  respond  to  toxicant  and  other  damaging  challenges 
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with  changes  in  transcription  up  to  10%,  making  them  an  ideal  model  for  this  study  5.  Further, 

C  JC  ?o 

this  model  has  already  been  used  for  both  sublethal  arsenic  and  low  dose  IR  studies  ’  ’  . 

Our  investigation  begins  to  address  the  gap  in  the  toxicology  literature  regarding  potential 
interactions  of  low  level  toxicant  exposures.  It  has  been  shown  that  transcriptional  changes 
within  the  cell  do  not  correlate  completely  with  translational  studies  '  .  Few  studies  have 
focused  on  the  proteomic  differences  induced  by  these  toxicants  and  this  study  begins  to  address 
the  lack  of  proteomic  data  and  provides  complementary  information  to  transcriptional  data.  As 
proteins  are  the  components  through  which  a  cell  enacts  change,  the  differences  found  in  the 
proteome  may  better  reflect  the  actual  cellular  response  of  significance  for  human  health 
outcomes.  Using  a  proteomics  approach  in  a  human  keratinocyte  model  to  mimic  human  skin 
exposure,  we  examined  the  interaction  of  ionizing  radiation  and  arsenic. 


Experimental  Procedures 

Cell  Culture 

A  spontaneously  immortalized  human  keratinocyte  (SIK)  cell  line  was  grown  with  a 
lethally  irradiated  feeder  layer  of  3T3  cells  obtained  from  ATCC  #48-X  (ATCC,  VA)  5.  Cells 
were  supported  with  a  3:1  mixture  of  Dulbecco-Vogt  Eagle  and  Ham’s  F-12  media  containing 
5%  fetal  bovine  serum,  0.4  pg/ml  hydrocortisone,  5  pg/ml  insulin,  5  pg/ml  transferrin,  20  pM 
triiodothyronine,  0.18  mM  adenine,  10  ng/ml  cholera  toxin  and  10  ng/ml  epidermal  growth 
factor 5.  Cell  medium  was  changed  every  3  days. 
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Arsenic  and  IR  Treatment 


Once  cells  reached  near  confluency,  half  of  the  flasks  were  treated  with  medium 
containing  2  pM  sodium  arsenite.  At  24  hours,  flasks  were  irradiated  with  x-rays  at  low  doses  of 
1  or  10  cGy.  Tissue  culture  flasks  were  irradiated  with  a  Varian  2100C  on  a  30  x  30  cm  acrylic 
block  with  a  tissue  equivalent  bolus  covering  the  top  of  the  flasks  to  ensure  accuracy  of  dose 
delivery.  Dose  rate  was  set  at  80  cGy/min,  and  SSD  was  101.3  cm  (1  cGy)  or  100  cm  (10  cGy). 
Radiation  free  controls  were  maintained  with  and  without  sodium  arsenite.  One  or  four  days 
post  irradiation,  flasks  were  rinsed  with  0.02%  EDTA  in  phosphate  buffered  saline  to  remove 
residual  3T3  cells.  All  samples  were  prepared  in  triplicate.  Flasks  were  held  at  -80°  C  until 
protein  extraction. 

Protein  Extraction 

Flasks  were  removed  from  -80°C  freezer  and  thawed  to  room  temperature.  One  ml  of 
Mammalian  Protein  Extraction  Reagent  (MPER,  Pierce  Biotechnology,  IL)  containing  1:100 
protease  inhibitor  (Sigma-Aldrich,  MO)  was  added  to  each  flask,  cells  were  scraped  and  lysate 
was  transferred  to  a  2  ml  tube  on  ice.  Each  sample  was  sonicated  for  1  min  and  centrifuged  at 
10,000  x  g  for  20  minutes  at  4°C.  Supernatant  was  removed  and  placed  in  a  clean  2  ml  tube. 
Protein  quantitation  was  performed  using  Coomassie  Plus  Protein  Assay  (Pierce  Biotechnology, 
IL). 

2D  electrophoresis 
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Isoelectric  focusing  was  performed  using  a  Protean  IEF  Cell  (Bio-Rad,  CA).  A  30  pg 
aliquot  of  the  protein  was  combined  with  a  lysis  buffer  solution  (0.5%  Triton  X-100,  4% 

CHAPS,  7  M  urea,  2  M  thiourea,  nanopure  water),  1%  Biolyte  3-10  buffer,  2%  protease  inhibitor 
cocktail  (Calbiochem,  CA),  0.065%  dithiothreitol  and  a  trace  amount  of  bromophenol  blue  dye 
for  a  total  volume  of  200  pi.  The  solubilized  protein  samples  were  left  at  room  temperature  for 
one  hour  before  loading.  ReadyStrip  IPG  strips  (pH  3-10, 1 1cm)  were  used  for  separation  in  the 
first  dimension  (Bio-Rad,  CA). 

Isoelectric  focusing  was  conducted  at  20°C  using  the  following  voltage  program:  50  V 
for  a  12-hour  rehydration  period;  50-250  V  linear  ramp;  250-8000  V  linear  ramp;  hold  (total  42 
kVh).  After  focusing,  the  strips  were  incubated  in  an  equilibration  buffer  (5  ml  consisting  of  50 
mM  Tris,  pH  8.8,  6  M  urea,  30%  glycerol,  2%  SDS,  trace  bromophenol  blue  and  0.065% 
ditheothreitol  (DTT))  for  1 5  minutes  on  a  rocking  platform.  The  strips  were  subsequently 
incubated  with  the  same  equilibration  buffer  substituting  10  mM  iodoacetamide  for  DTT  to 
alkylate  cysteine  sulfhydryls.  The  strips  were  then  placed  on  top  of  12%  SDS  Duracryl  gels  and 
sealed  using  0.5%  agarose  (Genomic  Solutions,  MI). 

The  second  dimension  separation  was  performed  in  a  Hoefer  SE  600/SE  660  2D-PAGE 
system  interfaced  to  a  Powerpac  3000  power  supply  (Bio-Rad,  CA).  Gels  were  run  in  an  SDS 
tank  buffer  (25mM  Tris,  192  mM  glycine,  0.1%  SDS)  at  15  mA  per  gel  for  30  minutes  followed 
by  25  mA  per  gel  until  the  dye  migrated  to  the  bottom  of  the  gel.  Broad  range  Precision  Plus 
Protein  Standard  plug  molecular  weight  protein  markers  were  used  for  mass  calibration  of  the 
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gels  (10-250  kDa)  (Bio-Rad,  CA).  Gels  were  fixed  in  10%  acetic  acid,  40%  methanol,  and  50% 
water,  silver  stained  and  scanned  with  an  Epson  Perfection  4870  photo  scanner. 

Image  analysis 

The  36  gel  images  were  processed  with  analysis  software  Progenesis  (PG240  v2006)  and 
TT900  S2S  (Nonlinear  Dynamics,  UK).  Gel  images  were  first  warped  with  TT900  S2S.  Warped 
images  were  then  imported  to  Progenesis  for  further  analysis  including:  spot  detection,  spot 
matching,  background  subtraction,  spot  filtering  and  ‘Samespot  Outline’.  The  Samespot  Outline 
in  Progenesis  copies  spot  outlines  from  gels  where  a  spot  exists  to  those  gels  missing  the  spot, 
then  calculates  the  spot  volume  within  the  new  outlines.  Therefore  all  missing  values  are  filled 
with  calculated  volumes.  There  were  2,002  individual  proteins  recognized  on  each  gel  with 
Progenesis  (PG240  v2006.  Nonlinear  Dynamics,  UK). 

Protein  digestion  and  mass  spectrometry 

The  Nevada  Proteomics  Center  analyzed  selected  proteins  by  trypsin  digestion  and 
MALDI  TOF/TOF  analysis.  Gel  pieces  were  destained  using  reagents  from  the  Silver  Quest  - 
Silver  staining  kit  (Invitrogen,  CA,  cat#  LC6070).  Spots  were  digested  using  a  previously 
described  protocol  with  some  modifications 33.  Samples  were  washed  twice  with  25  mM 
ammonium  bicarbonate  (ABC)  and  100%  acetonitrile,  reduced  and  alkylated  using  10  mM  DTT 
and  100  mM  iodoacetamide  and  incubated  with  75  ng  sequencing  grade  modified  porcine  trypsin 
(Promega,  WI)  in  25  mM  ABC  for  6  hours  at  37°C.  Samples  were  spotted  onto  a  MALDI  target 
with  ZipTipp-C18  (Millipore  Corp.,  MA).  Samples  were  eluted  with  70%  acetolnitrile,  0.2% 
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formic  acid  and  overlaid  with  0.5  pi  5  mg/ml  MALDI  matrix  (a-Cyano-4  hydroxycinnamic  acid, 
10  raM  ammonium  phosphate).  All  mass  spectrometric  data  were  collected  using  an  ABI  4700 
Proteomics  Analyzer  MALDI  TOF/TOF  mass  spectrometer  (Applied  Biosystems,  CA),  using 
their  4000  Series  Explorer  software  v.  3.6.  The  peptide  masses  were  acquired  in  reflectron 
positive  mode  (1-keV  accelerating  voltage)  from  a  mass  range  of  700  -  4000  Daltons  and  either 
1250  or  2500  laser  shots  were  averaged  for  each  mass  spectrum.  Each  sample  was  internally 
calibrated  on  trypsin’s  autolysis  peaks  842.51  and  221 1.10  to  within  20  ppm.  Any  sample  failing 
to  internally  calibrate  was  analyzed  under  default  plate  calibration  conditions  of  1 50  ppm.  Raw 
spectrum  filtering/peak  detection  settings  were  S/N  threshold  of  3,  and  cluster  area  S/N 
optimization  enabled  at  S/N  threshold  10,  baseline  subtraction  enabled  at  peak  width  50.  The 
eight  most  intense  ions  from  the  MS  analysis,  which  were  not  on  the  exclusion  list,  were 
subjected  to  tandem  mass  spectrometry.  The  MS/MS  exclusion  list  included  known  trypsin 
masses  along  with  unidentified  background  peaks:  842.51,  856.52,  870.54,  101 1.65,  1045.56, 
1126.56,  1338.83,  1666.01,  1794.9,  1940.94,  221 1.10,  2225.12,  2283.18  and  3094.62.  For 
MS/MS  analysis  the  mass  range  was  70  to  precursor  ion  with  a  precursor  window  resolution  of 
50  FWHM  (full-width  at  half  maximum)  with  an  average  2500  laser  shots  for  each  spectrum, 

CID  on,  metastable  suppressor  on.  Raw  spectrum  filtering/peak  detection  settings  were  S/N 
threshold  of  5,  and  cluster  area  S/N  optimization  enabled  at  S/N  threshold  6,  baseline  subtraction 
enabled  at  peak  width  50.  The  data  was  then  stored  in  an  Oracle  database  (Oracle  database 
schema  v.  3.19.0,  Data  version  3.90.0). 

MALDI  data  analysis 
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The  data  was  extracted  from  the  Oracle  database  and  a  peak  list  was  created  by  GPS 
Explorer  software  v  3.6  (Applied  Biosystems)  from  the  raw  data  generated  from  the  ABI  4700. 
Analyses  were  performed  as  combination  mass  spectrometry  and  tandem  mass  spectrometry. 

MS  peak  filtering  included  mass  range  700-4000  Da,  minimum  S/N  filter  10,  mass  exclusion 
tolerance  of  0.2  Da.  Exclusion  list  of  known  trypsin  fragments  and  unidentified  background 
peaks:  221 1.2,  2283.2,  1045.6,842.5,  1794.9,  1011.65,  1338.83,  1666.01.  A  peak  density  filter 
of  50  peaks  per  200  Da  with  a  maximum  number  of  peaks  set  to  65.  MS/MS  peak  filtering 
included  mass  range  of  60  Da  to  20  Da  below  each  precursor  mass.  Minimum  S/N  filter  5,  peak 
density  filter  of  50  peaks  per  200  Da,  cluster  area  filter  used  with  maximum  number  of  peaks  65. 
The  filtered  data  were  searched  by  Mascot  v  1 .9.05  (Matrix  Science)  using  CDS  combined 
database  (Celera  Discovery  System  v.  KBMS3. 2.200401 19),  containing  1,416,555  sequences. 
Searches  were  performed  without  restriction  to  protein  species,  Mr,  or  pi  and  with  variable 
oxidation  of  methionine  residues  and  carbamidomethylation  of  cysteines  (no  fixed 
modifications).  Maximum  missed  cleavage  was  set  to  1  and  limited  to  trypsin  cleavage  sites. 
Precursor  mass  tolerance  and  fragment  mass  tolerance  were  set  to  20  ppm  and  ±  0.2  Da, 
respectively.  Protein  hits  with  high  confidence  identifications  and  statistically  significant  search 
scores,  greater  than  95%  confidence  interval  (C.I.%)  or  p  <  0.05,  were  accepted.  High 
confidence  identifications  were  consistent  with  the  protein  experimental  Mr,  and  pi  and  the 
majority  of  the  ions  present  in  the  mass  spectra  were  accounted  for. 

Immunoblot  analysis 

Western  blots  were  performed  for  sequenced  proteins  that  had  high  confidence 
identification  and  which  had  available  commercial  antibodies.  Protein  samples  (10  pg)  were 
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separated  on  12%,  Ready  Gel  Tris-HCl  precast  gels  (BioRad,  CA).  Proteins  were  electroblotted 
to  PVDF  membrane  and  blocked  overnight  at  4°C.  Primary  antibodies  were  pyruvate  kinase 
(ab6 191,  Abeam,  MA),  a-enolase  (sc- 15343,  Santa  Cruz  Biotechnology,  CA),  S100A9  (sc-8114, 
SCB,  CA),  PDI  (sc-30932,  SCB,  CA),  profilin-1  (sc- 18346,  SCB,  CA),  annexin  XI  (sc-9322, 
SCB,  CA),  E-FABP  (sc- 16060,  SCB,  CA),  and  LDH-A  (sc-27230,  SCB,  CA),  cytokeratin  1  (sc- 
17091,  SCB,  CA),  CaM  (sc-1989,  SCB,  CA),  HSP27  (sc-1048,  SCB,  CA),  HINT-1  (10717-1- 
AP,  Proteintech  Group,  IL)  and  cyclophilin  A  (10720-1-AP,  Proteintech  Group  ,  IL).  All 
primary  antibodies  were  used  at  0.2  pg/ml  final  concentration,  typically  a  1:1000  dilution.  The 
secondary  antibody  (donkey  anti-goat-hrp,  sc-2020  or  donkey  anti-rabbit,  sc-2004,  SCB,  CA) 
was  used  at  a  1 :40,000  dilution.  Membranes  were  developed  in  ECL  Advanced  (GE  Healthcare, 
NJ)  and  images  were  captured  on  a  ChemiDoc  system  with  Quantity  One  software  (BioRad, 

CA).  A  monoclonal  antibody  to  P-actin  was  used  as  a  loading  control  and  all  density  readings 
were  normalized  (sc-47778,  SCB,  CA).  Western  blots  were  performed  in  triplicate.  The  western 
blot  data  were  analyzed  using  ANOVA  after  normalization  by  p-actin. 

Statistical  Analysis 

Statistical  differential  analyses  consisted  of  data  preprocessing,  analysis  of  variance 
(ANOVA),  adjusting  for  multiple  hypothesis  testing  and  the  chi-square  test.  Statistical  methods 
are  often  based  on  the  assumption  that  data  are  normally  distributed  with  constant  variance  not 
dependent  on  the  mean;  therefore,  it  is  important  to  preprocess  the  data  so  that  they  do  not 
violate  these  assumptions.  In  this  case  preprocessing  consisted  of  logarithmic  transformation  and 
mean-centering  normalization.  The  analysis  of  variance  model  was  applied  to  identify 
differentially  expressed  proteins  with  effects  for  the  level  of  arsenic,  the  level  of  IR,  and  the 
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interaction  effect.  Two  methods  of  estimating  the  protein-specific  variance  in  ANOVA  were 
utilized,  the  usual  mean  square  for  error  and  an  estimate  adjusted  by  an  empirical  Bayes  method 
originally  developed  for  microarrays  34, 35  (manuscript  in  progress,  Dan  Li). 

Results 

Image  and  Statistical  Analysis 

Proteins  were  isolated  from  keratinocytes  that  had  been  exposed  to  0  or  2  pM  sodium 
arsenite  and  0,  1  or  1 0  cGy  of  irradiation  for  one  or  four  days.  These  proteins  were  separated 
using  two-dimensional  gel  electrophoresis.  All  conditions  were  run  in  triplicate  and  the  resulting 
36  gels  were  imaged  and  analyzed.  The  initial  analysis  included  keratinocytes  exposed  to  a  1 
cGy  irradiation  dose.  Little  differential  expression  was  detected  with  this  low  dose  and  these 
samples  were  omitted  from  the  final  analysis. 

ANOVA  was  performed  on  the  each  of  2,002  spots  detected  at  each  time  point 
individually.  Proteins  that  displayed  significant  (p  <  0.05)  differential  expression  (after 
correction  for  multiple  comparisons)  using  either  ANOVA  method  were  selected  as  candidates 
for  sequencing,  resulting  in  444  spots  identified  for  further  characterization. 

Mass  Spectrometry  and  Protein  Identification 

Protein  spots  were  chosen  for  sequencing  only  if  they  were  distinct,  outside  of  areas 
with  background  smearing  and  could  be  removed  from  the  gel  without  excising  other  nearby 
proteins.  A  total  of  24  samples  were  sent  to  the  Nevada  Proteomics  Center  for  analysis,  nine 
spots  from  day  one  and  fifteen  from  day  four  (figure  1).  Thirteen  proteins  with  a  protein  score 
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CI%  above  95%  were  identified  by  mass  spectrometry  and  tandem  mass  spectrometry  (table  1 ). 
The  proteins  identified  by  mass  spectrometry  and  tandem  mass  spectrometry  from  the  24  hour 
time  point  included:  calmodulin  (CaM),  heat  shock  protein  27  (HSP27),  lactate  dehydrogenase  A 
(LDH-A)  and  protein  disulfide  isomerase  precursor  (PDI,  synonym:  thyroid  hormone  binding 
protein  precursor)  (figure  2).  The  proteins  found  four  days  post  exposure  included:  annexin  XI 
(annexin  XI),  S100A9  (synonym:  calgranulin  B),  cyclophilin  A  (synonym:  peptidyl-prolyl  cis- 
trans  isomerase),  a-enolase,  epidermal  fatty  acid  binding  protein  (E-FABP),  histidine  triad 
nucleotide-binding  protein  1  (HINT-1),  profilin-1,  pyruvate  kinase  M  isozyme,  and  R3372_l  (a 
protein  fragment)  (figure  3). 

Immunoblotting 

Twelve  identified  proteins  had  available  commercial  antibodies  and  were  selected  for 
immunoblotting  (table  2).  This  was  done  as  an  important  cross-check  on  the  ANOVA  analysis 
for  the  2D  gels.  Given  the  complexity  of  the  required  statistical  analysis  for  the  2D  gels,  it  is 
always  possible  that  some  spots  identified  as  differentially  expressed  were  in  fact  artifactual. 
CaM  was  not  recognized  by  the  antisera  in  detectable  amounts  and  no  protein  band  was  verified 
at  the  appropriate  molecular  weight.  This  does  not  exclude  the  possibility  that  CaM  may  have 
been  differentially  expressed  since  western  blots  require  a  relatively  high  concentration  for 
detection.  For  the  eleven  remaining  proteins,  the  P-actin-normalized  intensity  data  determined 
(table  3).  These  data  were  analyzed  using  ANOVA  (table  4).  Eight  of  the  eleven  proteins  (E- 
FABP,  a-enolase,  HINT-1,  HSP27,  LDH-A,  PDI,  cyclophilin  A,  and  S100A9)  showed 
significant  (p  <  0.05)  differences  for  either  IR,  As,  or  the  interaction,  which  far  exceeds  the 
chance  rate  (roughly  one  or  two  false  significance  values  could  be  expected  out  of  the  1 1 
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proteins  by  chance  alone).  Three  showed  no  significant  difference  at  the  p  <  0.05  level  (annexin 
XI,  pyruvate  kinase  and  profilin-1).  The  HSP27  was  down  regulated  and  HINT-1  was  up 
regulated  in  response  to  individual  treatments,  though  neither  protein  showed  an  interaction 
effect  that  was  significantly  different  than  either  single  treatment.  Cyclophilin  A,  E-FABP,  PDI 
and  S100A9  decreased  with  exposure  to  arsenic.  Both  E-FABP  and  LDH-A  showed  a  response 
to  the  combined  insult  that  was  greater  than  would  have  been  expected  from  either  treatment 
alone,  while  a-enolase  and  PDI  had  a  response  that  was  less  than  would  be  expected  from  the 
individual  exposures. 

Discussion 

Low  dose  ionizing  radiation  and  low  levels  of  arsenic  are  environmental  toxins  and 
below  a  certain  threshold,  their  individual  toxicity  is  believed  to  be  minimal.  However  the 
effects  of  combining  these  two  known  carcinogens  at  these  doses  are  unknown.  As  a  model  for 
low  dose  toxicant  interaction,  they  are  of  significant  interest  as  irradiation  is  ubiquitous  and 
arsenic,  although  heavily  regulated  in  the  US,  is  still  a  major  environmental  burden  in  many 
countries.  There  is  little  in  the  scientific  literature  examining  human  health  risks  associated  with 
possible  toxicant  interactions  despite  the  known  synergy  of  such  agents  in  malignancy  ',2’ l5'36. 
Thus,  these  data  serve  to  begin  to  fill  this  information  gap. 

Of  the  twelve  proteins  selected  for  immunoblotting,  eight  showed  significant  differences 
in  protein  expression  and  three  yielded  inconsistent  results  (annexin  XI,  pyruvate  kinase  and 
profilin-1).  It  is  possible  that  the  expression  changes  detected  by  the  ANOVA  of  the  gel  images 
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were  too  small  to  be  detected  by  western  blotting  given  its  semi-quantitative  nature.  All  proteins 
with  variable  results  had  faint  spots  on  the  2D  gels  indicating  a  very  low  protein  concentration  as 
silver  staining  can  detect  nanogram  quantities  of  protein,  a  sensitivity  beyond  that  of  western 
blotting. 

Proteins  that  were  immunoblotted  included  three  proteins  from  the  early  time  point:  PDI, 
LDH-A  and  HSP27.  None  of  these  proteins  were  found  at  the  later  time  point  indicating  a 
transient  response.  At  four  days  post  exposure,  E-FABP,  a-enolase,  S100A9,  HINT-1  and 
cyclophilin  A  had  altered  expression.  Doses  examined  in  this  study  are  much  lower  than  those 
that  are  clearly  cytotoxic  and  this  proteomic  response  implies  that  the  cells  are  actively 
responding  to  low-level  exposures. 

Our  data  adds  to  and  confirms  several  prior  studies  using  these  toxicants.  PDI,  a  stress 
response  protein  was  down  regulated  in  a  previous  study  by  exposure  to  arsenic  5.  We  found  this 
protein  was  down  regulated  in  response  to  both  arsenic  and  arsenic  plus  IR.  The  LDH 
isoenzyme  spectrum  in  human  serum  changes  following  radiation  exposure,  tumor  development, 
as  well  as  in  response  to  disease  37'39.  The  isoform  LDH-A  is  typically  found  to  be  down 
regulated  in  response  to  these  stressors.  This  protein  was  identified  as  significantly  down 
regulated  in  response  to  the  combined  exposure.  However,  the  western  blot  images  and  the 
normalized  density  readings  show  that  the  LDH-A  appears  to  be  down  regulated  in  response  to 
both  treatments  individually,  but  these  results  were  not  statistically  significant. 
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This  study  revealed  that  the  interaction  effect  of  IR  with  arsenic  significantly  down 
regulated  a-enolase.  Enolase  is  multifunctional  and  can  serve  as  a  plasminogen-binding 
receptor,  a  stress  protein,  a  structural  protein,  a  glycolytic  enzyme  and  a  transcriptional  repressor 
40,41.  Wu  et  al  reported  that  cells  treated  with  all-trans  retinoic  acid  (a  suppressor  of  cellular 
growth  and  differentiation)  showed  a  two-fold  decrease  in  a  -enolase  levels  in  hepatocarcinoma 
cells  and  suggested  that  a-enolase  is  down  regulated  in  slow  growth  conditions 42. 

Interestingly,  a  study  by  Tallman  et  al  shows  that  the  use  of  all-trans  retinoic  acid  in  the 
treatment  of  patients  with  acute  promyelocytic  leukemia  results  in  increased  survival  43 .  a- 
enolase  is  known  to  have  higher  expression  in  many  tumor  types  and  has  been  shown  to  be 
hypoxia  inducible,  however  caution  must  be  used  when  extrapolating  results  from  cancerous 
tissues  and  cell  lines  to  normal  human  keratinocytes  44' 45 .  It  has  been  postulated  that  a-enolase 
is  critical  to  tumor  invasion  and  metastasis  via  activation  of  the  plasminogen  system  41 .  While 
no  firm  conclusions  can  be  drawn,  this  interaction  of  low  dose  IR  and  low  dose  arsenic  in  a- 
enolase  expression  does  suggest  that  further  study  is  warranted  given  the  postulated  role  in 
malignant  behavior  and  in  light  of  current  cancer  therapies. 

Expression  of  the  S 1 00 A9  protein  was  significantly  down  regulated  by  the  arsenic  or  the 
combination  treatment.  S100A9  is  one  of  a  family  of  SI  00  proteins  found  in  a  wide  variety  of 
cells  and  involved  in  calcium  signal  transduction,  cell  cycle  progression  and  cellular 
differentiation,  especially  in  keratinocytes  and  epidermal  tissues 46.  Two  microarray  analyses  on 
the  effects  of  low  dose  ionizing  radiation  (0-100  cGy)  on  full  thickness  human  skin  show  that,  as 
a  large  family,  the  SI 00  gene  group  is  down  regulated  within  the  first  24  hours  5’ 26' 47 .  Another 
recent  proteomic  analysis  of  laryngeal  cancer  found  decreased  expression  in  the  S100A9  protein 
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as  compared  to  nearby  healthy  tissue  48.  To  our  knowledge,  there  are  no  reports  on  S100A9 
arsenic  altered  expression. 

E-FABP  was  significantly  down  regulated  by  arsenic  and  up  regulated  by  the  combined 
exposure.  E-FABP  has  an  affinity  for  binding  to  fatty  acids  and  is  involved  in  keratinocyte 
differentiation,  synthesis  of  phospholipids  and  has  a  role  as  an  antioxidant  which  scavenges 
reactive  lipids  49, 50.  E-FABP  is  known  to  interact  with  S100A7  and  these  proteins  are  up 
regulated  in  psoriasis  and  certain  cancers  ’  .  In  prior  genomic  studies,  S100A9  and  E-FABP 
have  been  shown  to  be  up  regulated  in  response  to  tissue  damage  53’ 54.  While  these  interactions 
need  to  be  clarified,  our  data  adds  supports  the  possibility  that  these  two  proteins,  S100A9  and  E- 
FABP,  may  interact  with  one  another. 

Two  proteins  yielded  variable  immunoblotting  results,  pyruvate  kinase  M1/M2  isozyme 
and  cyclophilin  A.  The  polyclonal  pyruvate  kinase  antibody  was  not  specific  for  the  M 
isoform.  A  shift  in  the  pyruvate  kinase  isozyme  profile  from  tissue  specific  pyruvate  kinase  to 
the  M2  form  is  an  indication  of  early  carcinogenesis  and  the  conformational  shifting  of  M2  from 
tetrameric  to  dimeric  isoform  occurs  in  late  carcinogenesis  55'57.  The  ANOVA  analysis  of 
cyclophilin  A  revealed  a  significant  response  with  exposure  to  arsenic.  Cyclophilin  A  is 
involved  in  Ca2+  ion  signaling  and  when  complexed  with  cyclosporin  A  it  inhibits  calcineurin,  a 
calcium/calmodulin  phosphatase  58.  We  identified  calmodulin  and  PDI  regulation  upon 
exposure  to  arsenic,  while  a  prior  study  of  the  response  of  keratinocytes  to  arsenic  identified 
calcineurin,  all  of  which  are  involved  in  Ca  ion  signaling  ’  . 
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Both  IR  and  arsenic  dose  responses  are  currently  based  on  the  idea  of  a  linear  no¬ 
threshold  model  13,24.  Many  studies  are  revealing  that  low  dose  exposures  are  fundamentally 
different  from  high  dose  exposures  ’  .  For  example,  it  has  been  shown  that  a-enolase  was  up 
regulated  in-vivo  with  a  high  dose  exposure  of  9  Gy,  whereas  our  study  did  not  confirm  this 
trend  59.  Prasad  et  al  showed  that  a  dose  of  6  Gy  of  irradiation  led  to  increases  in  the  protein 
levels  of  PDI,  calreticulin  and  calnexin  in  apoptotic  cells  48  hours  post  exposure.  Low  dose 
exposure  of  2  pM  arsenic  alone  or  with  10  cGy  IR  diminished  cellular  amounts  of  both  PDI  and 
calreticulin  5.  These  results  suggest  that  the  response  to  low  doses  of  irradiation  and  arsenic  may 
be  substantially  different  than  those  seen  following  higher  dose  exposures.  A  response  pattern  of 
that  type  would  be  consistent  with  emerging  data  from  other  studies  dealing  with  low  dose  IR 
and  low  dose  arsenic  5’ 26, 47. 

Several  of  the  identified  proteins  were  found  to  have  a  response  to  the  combined 
exposure  that  was  different  from  the  response  to  either  toxicant  individually,  including  E-FABP, 
a-enolase,  LDH-A  and  PDI.  These  data  demonstrate  a  possible  response  to  the  combined  insult 
that  is  significantly  different  than  either  treatment  alone,  making  them  candidates  for  further 
study  as  potential  biomarkers.  Many  of  the  proteins  found  in  this  study  are  currently  recognized 
as  biomarkers  of  disease  processes.  Pyruvate  kinase  M2  is  a  biomarker  for  colorectal  cancer, 
while  a  change  in  the  usual  LDH  isozyme  spectrum  is  indicative  of  ischemia,  radiation  treatment 
or  cancer,  while  S100A9  has  been  found  as  a  serum  component  after  radiation  38,60,61.  The 
identification  of  biomarkers  in  conjunction  with  changing  isozyme  spectra  that  is  unique  to  a 
combination  of  environmental  toxicants  may  lead  to  novel  detection  panels  for  suspected 
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environmental  toxicant  exposures.  With  these  findings,  we  have  begun  the  process  of  identifying 
potential  biomarkers  of  these  types  of  combined  exposures. 

This  study  examined  near  background  levels  of  arsenic  and  IR  to  determine  the  potential 
for  overlapping  cellular  responses.  Previous  work  has  shown  that  arsenic  exposure  alters  the 
expression  of  heat  shock  proteins,  annexin  and  PDI,  while  IR  has  been  reported  to  alter  heat 
shock  protein  60,  annexin  and  PDI,  as  well  as  a-enolase  4’ 6’ 59.  Our  research  indicates  that 
HINT-1  and  HSP27  respond  in  the  same  manner  to  either  arsenic  or  IR.  Other  proteins  (E- 
FABP,  a-enolase,  LDH-A  and  PDI)  had  an  effect  that  was  significantly  different  than  would  be 
expected  based  on  a  single  exposure.  These  proteomic  responses  are  consistent  with  the  idea  that 
cellular  defenses  to  environmental  toxicants  have  partially  overlapping  response  pathways, 
suggesting  the  possibility  that  the  alteration  of  two  minimally  toxic,  sublethal  exposures  may 
combine  into  a  lethal  or  carcinogenic  response. 
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Figure  Legends 


Table  1 

Protein  identification.  Proteins  identified  from  mass  spectrometry  and  tandem  mass 
spectrometry.  Fields  include:  spot  number  as  displayed  in  figure  1 ,  protein  name,  accession 
number  from  CDScombined  database  (Celera  Discovery  System  v.  KBMS3. 2.200401 19), 
MASCOT  protein  score,  MASCOT  ion  score,  protein  and  ion  score  confidence  indices,  species 
of  protein  hit,  and  protein  molecular  weight,  pi,  peptide  count  and  total  amino  acid  sequence 
coverage.  Also  shown  are  the  matched  peptides  and  corresponding  masses,  mass  errors, 
sequences  and  number  of  unmatched  sequences  if  protein  ion  score  is  less  than  100.  Peptides 
analyzed  by  tandem  mass  spectrometry  display  ion  scores  and  confidence  indices. 

Table  2 

Immunoblot  images  of  proteins  identified  by  MS/MS.  Immunoblots  were  performed  for  each 
selected  protein.  1 0  pg  of  protein  were  separated,  transferred  to  PVDF,  and  developed  using 
ECL  Advanced  reagent.  Images  were  captured  on  a  ChemiDoc  system.  (A)  control  sample  (B) 
IR  only,  (C)  As  only  and  (D)  IR  plus  As.  Each  immunoblot  was  performed  in  triplicate  and 
images  shown  are  representative  of  all  blots  performed.  Each  blot  was  probed  for  P-actin  as  a 
loading  control  and  one  representative  image  is  shown. 

Table  3 

Western  blot  density  data.  The  average  intensity  per  mm2  was  determined  for  each  band  using 
Quantity  One  software.  The  intensity  data  from  each  protein  were  normalized  to  the  P-actin 
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data.  Cyclophilin  A  yielded  variable  results  with  immunoblotting,  but  was  shown  to  have  a 
significant  response  to  arsenic  in  the  ANOVA  analysis.  The  data  from  three  replicates  blots  are 
given. 

Table  4 

Significant  effects  in  the  ANOVA  of  the  western  blot  analysis.  The  p-value  is  given  for  all 
effects  significant  at  the  0.05  level.  The  (T)  symbol  in  the  IR  or  As  column  indicates  that  the 
given  protein  is  up  regulated  upon  treatment  with  the  respective  toxicant,  and  the  (i)  symbol 
indicates  down  regulation.  For  the  IR*As  interaction,  the  (+)  symbol  indicates  that  the  protein 
expression  is  greater  than  would  have  been  expected  from  the  two  stimuli  independently,  and  (-) 
indicates  that  the  protein  expression  is  less  than  would  have  been  expected.  In  general,  we  can 
conclude  that  there  has  been  evidence  of  an  IR  perturbation  of  protein  expression  if  either  the  IR 
or  the  IR*As  effect  is  significant,  and  that  there  has  been  evidence  of  an  As  perturbation  if  either 
the  As  or  IR*As  effect  is  significant.  Thus  6  of  the  1 1  proteins  showed  a  reaction  to  IR,  8  of  the 
1 1  showed  a  reaction  to  As  (including  all  that  showed  a  reaction  to  IR),  and  3  were  inconclusive. 
The  ANOVA  was  performed  for  all  exposures,  —  indicates  no  significant  change. 

Figure  1 

Figure  1  Two  dimensional  gels  with  labeled  protein  spot  numbers  selected  for  analysis. 

Keratinocytes  were  treated  with  low  dose  arsenic  and  ionizing  radiation  and  protein  was  isolated 
one  or  four  days  post  exposure.  Numbered  spots  were  excised  and  sent  to  the  Nevada  Proteomic 
Center  for  identification. 
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Figure  2 

MS  spectrum  of  proteins  identified  at  the  24  hours  time  point.  The  images  are  annotated  in 
color  to  correspond  with  the  MASCOT  results.  (A)  CaM,  (B)  HSP27,  (C)  LDH-A  and  (D)  PDI. 
Red  annotation  denotes  a  direct  match  to  the  protein  in  MASCOT,  blue  indicates  submission  to 
MASCOT  with  no  match.  Masses  shown  in  black  were  excluded  from  the  MASCOT  search  as 
they  were  a  trypsin  peak,  a  known  contaminant  or  the  S/N  was  lower  than  the  threshold. 

Figure  3 

MS  spectrum  of  proteins  identified  fours  days  post  IR  treatment.  The  images  are  annotated 
in  color  to  correspond  with  the  MASCOT  results.  (A)  annexin  XI,  (B)  cyclophilin  A,  (C)  a- 
enolase,  (D)  HINT-1,  (E)  E-FABP  (spot  2611),  (F)  E-FABP  (spot  2640),  (G)  profilin-1,  (H) 
pyruvate  kinase  M  isozyme,  (I)  R3372_l  and  (J)  S100A9.  Red  annotation  is  a  match  to  the 
protein  hit  in  MASCOT,  blue  labeling  indicate  submission  to  MASCOT  with  no  match.  Black 
indicates  masses  excluded  from  the  MASCOT  search  as  they  were  trypsin  peaks,  known 
contaminants  or  the  S/N  was  lower  than  the  threshold. 
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